Hartree-Fock ground state of the two-dimensional electron gas with Rashba spin-orbit 

interaction 
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We search for the uniform Hartree-Fock ground state of the two-dimensional electron gas formed 
in semiconductor heterostructures including the Rashba spin-orbit interaction. We identify two 
competing quantum phases: a ferromagnetic one with partial spin polarization in the perpendicular 
direction and a paramagnetic one with in-plane spin. We present a phase diagram in terms of 
the relative strengths of the Rashba to the Coulomb interaction and the electron density. We 
compare our theoretical description with existing experimental results obtained in GaAs-AlGaAs 
heterostructures. 

PACS numbers: 71.10.Ca, 71.70.Ej, 73.21.-b 



I. INTRODUCTION 

The two-dimensional electron gas (2DEG) is a paradig- 
matic system of semiconductor physics and technology. 
Traditionally, the electronic spin degree of freedom in 
this system has played a secondary role. This situa- 
tion has changed recently with the emergence of the 
promising field of spintronics.^'^ The most controllable 
and often predominant spin-orbit coupling in semicon- 
ductor 2DEGs is the Rashba interactioni^ It is then 
important to determine the various many-body prop- 
erties of the 2DEG in its presence. The ground state 
of the uniform 2DEG without spin-orbit interaction is 
not known exactly^-®'^ In this Brief Report we con- 
centrate on how the single-particle Rashba term affects 
the ground state in the Hartree-Fock (HE) mean-field 
approximation.® While for the uniform electron gas with- 
out spin-dependent potentials, the HE approach yields 
trivial single-particle spin-orbitals (plane waves and pure 
spin states), the presence of the Rashba spin-orbit in- 
teraction causes the HE solution to possess an intriguing 
spin texture in momentum space. In this Brief Report we 
formulate the HE theory and obtain solutions for the uni- 
form 2DEG with Rashba interaction. Our main finding 
is a spatially uniform ferromagnetic phase, characterized 
by a net out-of-plane partial magnetization appearing in 
a window of densitiesi^ 



II. HARTREE-FOCK THEORY WITH RASHBA 
SPIN-ORBIT INTERACTION 

The Hamiltonian of the 2DEG in the presence of the 
Rashba spin-orbit interaction is 



H 



E^«.^ + |E-(^ 



(1) 



«#i 



The second term gives the interparticle Coulomb inter- 



action v{ri — Tj 



/e|ri — Tjl, where e is the dielectric 



sum of the Rashba Hamiltonians of the individual elec- 
trons. The latter are given by 



2m* \ dy ox 



(2) 



Here, m* is the conduction-band effective mass, a{Ez) is 
a structural parameter that determines the strength of 
the Rashba coupling, and Ux and Uy are Pauli matrices. 
The Rashba Hamiltonian can be solved analytically, and 
the following wave functions and energies are obtained)^ 



V'ks(r) = 



E{k,s) 



1 



2A 



2m* 



sa{Ez)k. 



(3) 
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In these expressions k is the two-dimensional wave vec- 
tor, if is the angle of k in polar coordinates, and A is the 
surface area of the sample. The spin quantum number 
s = ±1 denotes spin-up and spin-down eigenstates with 
respect to the spin quantization axis, which lies in the 
x-y plane and is perpendicular to k with a polar angle 
4>R(k) = If — 7r/2. Notice that, however, the Rashba 
ground state is paramagnetic)^ Also, the Rashba Hamil- 
tonian Hr is time-reversal invariant. This invariance 
requires that V'ks(r) and '0-ks(r) be Kramers-conjugate 
states with the same energy eigenvalue Eik^ s) (Ref. llOf ). 

Eor a spatially uniform solution the HE spin-orbitals 
can be written as 



V'ks(r) = 



1 



Mk,s) 



(5) 



constant of the semiconductor, and the first term is the 



The spinor x(k, s) has components x±0^is)i which 
arc the unknown spinor amplitudes to be determined 
through the HE procedure. Notice that while the Rashba 
amplitudes of Eq. ([3]) depend only on the polar angle f 
of k, we allow x±(k, s) to depend also on the modulus of 
k. The spin quantum number s = ±1 denotes, like in the 



noninteracting Rashba problem, the up- and down-spin 
eigenstates in an unknown spin-quantization axis u(k) 
with polar angles 9(k) and 0(k). 

The functional that one has to minimize is given by (di- 
rect Coulomb terms drop out in the jellium-model elec- 
tron gas) 

T = gx^(k„ s)[i?fl(k,) - i?(k„ s)]x(k„ s) 

s,i 



2A 



^ 7;(k,-k,)[xUk„s')x(k.,s)] 



x[x^(kj,s)x(kj,s')], 



(6) 



where Ns are the numbers of occupied orbitals with s = 
±1, £'(ki,s) are the HF single-particle energies, w(ki — 
kj) = 27re^/£|ki — kj| is the Fourier transform of the 
Coulomb interaction, and 



/ h^P/2m* ia{E^)ke-'^v \ 



HR{k) 



(7) 
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We minimize this functional with respect to the ampli- 
tudes x*(ki, s) and obtain the single-particle energies, 

E{k, s) = h sa(^^)fc sin6'(k) 



2m* 



1 
'2A 



J2 t;(k-k')[l + ss'u(k)-u(k')]. (8) 



s'.k'GP, 



By demanding that the single-particle HF Hamilto- 
nian be diagonal on the basis of the spinor amplitudes 
Xe(ki,±) we obtain the following integral equations: 



2a{E^)ksm[(l){k) - ^]cos6l(k) = 
J Y, s't;(k-k'){sin6'(k)cos6l(k') 

s'.k'Gl',/ 

-sin6l(k')cos6'(k)cos[</>(k) 



A 



(k')]}, 



(9) 



2a{E^)kcos[(j>{]<i)-ip]^^ Y^ s'u(k-k') 

X sin6i(k') sin[<?!>(k) - </.(k')]. (10) 

In these equations, we have substituted the following ex- 
pressions for X£(ki,±) in order to display the depen- 
dence on the polar angles 6'(k) and ^(k): x+(k, +) — 
cos(6'/2)cxp(-i(^/2), x-(k,+) = sin(6'/2) exp(i(/)/2), 
X+(k,-) = -sin(0/2)exp(-z(^/2), and X-(k,-) - 
cos{9 /2) exp(i0/2). The summation domains 'D± are the 
regions of k space occupied by electrons and their ar- 
eas are N±, respectively. Notice that Eq. ([5]) reduces to: 
(i) Eq. (dl when the Coulomb interaction is neglected 
because the third term of the right-hand side (RHS) of 
Eq. ([5]) drops out and 6{k) = 7r/2 in the single-particle 
Rashba problem; (ii) the HF one-particle energy if the 
Rashba coupling is omitted because a = and u(k) = z. 



III. ISOTROPIC SOLUTION 

In the absence of Pomeranchuk instabilitieaii' (PI) (de- 
formations of the Fermi sphere in three dimensions or cir- 
cle in two dimensions), the domains 'D± should be taken 
as having circular symmetry. The issue of the occur- 
rence of PI in Fermi liquids with isotropic central interac- 
tions is currently being studied in two-dimensional (2D) 
and three-dimensional (3D) systemsi ^^'^'^ In the partic- 
ular case of the bare Coulomb interaction (which is our 
case) , the existing theory is not able to categorically pre- 
dict or rule out the PI. However, a screened Coulomb 
interaction does not produce PI (Ref. |13I), which can be 
taken as an indication that no PI occurs either for the 
bare interaction. In general, it is safe to assume that 
no PI will occur unless: (i) there is a well-defined length 
scale in the interaction, and (ii) that length scale is larger 
than the mean interparticle distancciii^ These conditions 
are clearly not satisfied by the bare Coulomb interac- 
tion, which does not possess a characteristic length scale. 
Based on these results, it is safe to assume that our inte- 
gration domains have circular symmetry. 

The circular symmetry of the integration domains im- 
plies that the dispersion relations E(k, s) in Eq. ^ must 
be isotropic. This in turn, requires that 0(k) be indepen- 
dent of if and that (/>(k) — ip — ix 12 as in the noninter- 
acting Rashba problem»i^ Thus Eq. (fTU]) is automatically 
satisfied and Eq. ^ becomes 



x' dx' dip' 



pxcos6{x) 

'xa Jq \I X^ + 2:'2 — 2xx' C0S(93 — (y9') 

x[sin6'(a;)cos0(a;') — cos6'(a;) sin0(x') cos((^ — (/3')].(11) 

We have introduced the parameter p — 2a{Ez)e/e^ 
which indicates the relative strength of the Rashba and 
Coulomb interactions. The integration limit x^ con- 
tains the information of the integration domains 'D± 
introduced earlier. Having in mind the dispersion re- 
lation of the noninteracting Rashba problem described 
in Eq. ([4]), the situation can be summarized as follows: 
when both branches are occupied (high density) , the do- 
mains 'D± are filled Fermi circles of radii kp± related by 
AttUs = fc|.^ + fc|,_, with ns = {N+ + N^)/A. In this 
case Xc — kp+fkp- and, in principle, it is free to vary 
from zero to one. If only the lower branch is occupied 
(low density), there are two possibilities. If the lower 
branch has a minimum at A: = (unlike the noninter- 
acting Rashba problem), then there is a gap at fc = 
between the two branches. In this case we define Xc 
as before, obtaining Xc = 0. The formation of a gap 
is enabled by the lifting of Kramers degeneracy at the 
single-particle level, due to the appearance of a sponta- 
neous magnetization. Kramers' theorem states that the 
degeneracy of half-integer spins can only be removed by 
a magnetic field; we return to this point in the discussion 
of Eq. (J12p . If the lower branch does have a minimum 



at fC ^ U we take X(, — '^min/kmaxi where Krnin (kmax) 

are the inner (outer) Fermi radii of the hollow circular 



domain, now being inus = k^ax ~ ^min- I^ this case, 
there is no gap at fc = 0, no (ground-state) magnetic mo- 
ment appears and the time-reversal symmetry (Kramers 
degeneracy) is preserved. In Eq. pTj) . a; = |x| = |k|//ci?_ 
or a; = |k|/A:mQ2;, according to the context. 

Clearly, the noninteracting Rashba states with 6 = 
7r/2, are a solution of Eq. (|11|) . We call this the in-plane 
(IP) paramagnetic phase, which is the one we mentioned 
as having the lower dispersion-relation branch with a 
minimum at fc ^ 0. The nontrivial solution with varying 
0{x) gives rise to an OP ferromagnetic phase where, in 
turn, the lower branch has its minimum at fc = 0. We 
solve Eq. (fTTj) starting with the initial guess of 9o{x) — 0, 
which gives Oi{x) after integration. With di{x) as an in- 
put, we obtain 02{x) and so on. We consider that con- 
vergence is achieved when 9n{x) and 9n~i{x) differ in 
less than 0.1 percent. Let us enumerate our main find- 
ings on 9{x): (i) 9(0) ~ for all values of p (the Rashba 
spin-orbit interaction causes no effect when k — and 
consequently, the spin quantization axis must lie in the 
z direction), (ii) 9{x) is a monotonically increasing func- 
tion for all values of p, (iii) 9{x) increases with increasing 
p, and (iv) 9{x) never crosses the value n/2. 

The gap at fc = between the two branches may be 
obtained from Eq. ([8]) by means of the above-mentioned 
properties, and it is given by 



AE^E{0,+)-E{0,-) 



cos9{x)dx. (12) 



In the IP phase 9{x) = 7r/2 and then A^ = 0. In the 
OP phase, in turn, we have AE ^ 0. This gap, far from 
being a peculiarity of the OP phase, appears also in the 
HF theory of the 2DEG without Rashba coupling. In 
fact, setting 9{x) = we get AE — {e^ /e){kF- — kF+) 
and the gap is unfailingly related to a polarized ground 
state, where kF+ ^ kp- (i.e., N+ ^ 7V_) (an identical 
expression for AE multiplied by 2/7r can be obtained in 
a 3DEG: see the expression for e^ when fc ^ on p. 82 
of Ref.i). 



IV. HARTREE-FOCK GROUND STATE 



The remaining task is, for given p and r^ (r^ = 
l/agy'TTOs, where a*^ = h^e/'m*e'^ is the effective Bohr 
radius), to determine which of these two phases has lower 
energy. For each phase and given r^ , the value of Xc that 
minimizes the energy is found numerically. In Fig. [1] 
we show a phase diagram in terms of the parameters 
rg and p. The striking feature of this diagram is that 
for given < p ^ 1.3, the OP phase appears within 
a window of densities. As expected, for p = we re- 
cover the paramagnetic-ferromagnetic transition of the 
2DEG HF approximation at r^ = 2.01 (Refs.[lland[i3). 
As p increases, the left transition moves slightly towards 
smaller rg; in other words, the presence of the Rashba 
coupling favors a spin-polarized phase, albeit this polar- 
ization is partial for nonzero p (the original HF ferro- 
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FIG. 1: Ground-state phase diagram in the Hartree-Fock ap- 
proximation in terms of the density parameter Vs and the 
Rashba to Coulomb energy ratio p. 




FIG. 2: Mean value of the perpendicular spin projection per 
particle {fi = 1) as a function of the Rashba to Coulomb 
energy ratio in the OP phase. 



magnetic phase has full polarization). Also the window 
of densities in the OP phase shrinks as p increases. The 
right transition originates in the fact that the system 
diminishes the ground-state energy by filling the lower 
branch that has minimum at fc 7^ (IP phase), when r^ 
increases at fixed p. 

The partial spin polarization of the (ferromagnetic) OP 
phase can be seen in Fig. [2l where we plot the mean value 
of S^ per particle ((5^) = {Sy) = for both IP and OP 
phases) given by (?i = 1), 



(5.) = 



xcos9{x)dx. 



(13) 



Xa=0 



Notice that {Sz) does not depend on r^ since Xc = for 
the OP phased 



V. COMPARISON WITH EXPERIMENT AND 
CONCLUSION 

In a recent experimental work, Ghosh et alM' report a 
possible spontaneous spin polarization in mesoscopic two- 
dimensional systems. They studied 2DEGs in asymmet- 
ric Si (5-doped GaAs/AlGaAs heterostructures with den- 



sities as low as n^ = 5 x lO^cm"^ (rg = 7.6). The temper- 
ature was set at T « 40 mK or equivalently T/Tp Ri 0.02 
since Tp = 2.3 K at r^ = 7.6. According to their interpre- 
tation of the data, these authors found partial spin polar- 
ization estimated as C = {N+ — N-)/N « 0.2 appearing 
in a window of densities of width Ar^ « 1.8, centered 
around r^ « 6.5. After a detailed analysis of their data, 
the authors rule out the Rashba coupling and invoke an 
exchange-driven spontaneous spin polarization in order 
to explain their observed split zero-bias peak (ZBP) in 
the differential conductance. On the other hand, they 
explain the fact that the spin polarization is partial as 
a finite-temperature effect. Recall that at zero temper- 
ature, this transition is of first order and is predicted 
to occur between Vs « 13 (Ref. '^ and rg ~ 25 (Rcf. 
0). However, this contradicts the theoretical finding of 
Dharma-wardana and Perrot,^° in the sense that partial 
spin polarization due to the exchange happens for T /Tp 
between 0.3 and 1.6, i.e., well above T/Tp « 0.02 as 
reported by Ghosh et al. 

Our theory supports the interpretation of the experi- 
ment of Ghosh et al. based on an exchange-driven spin- 
polarized ground state. On the other hand, we explain 
the partial spin polarization in terms of the Rashba inter- 
action rather than in terms of finite temperature. This 
hypothesis could be tested by additional experiments 
with a tunable Rashba coupling, which would allow vary- 
ing our parameter p. In such experiments the degree of 
spin polarization could be modified and measured. 

Ghosh et al. observe a split ZBP within a window of 
electron densities. At the lower end of this window, 
disorder-induced localization destroys the electron-gas 



scenario. In a defect-free sample, we also predict a disap- 
pearance of the split ZBP but due entirely to the compe- 
tition between Rashba and exchange interactionsi^ At 
the higher end of the density window, the reason for the 
disappearance of the split ZBP is not as clear, and it is 
not specified by the authors. According to our theory, 
this feature could be interpreted as the ferromagnetic- 
paramagnetic transition of the electron gas at high den- 
sity, seen in Fig. [T] for low Vg. (Of course, the critical 
value of Vg is overestimated by the Hartree-Fock approx- 
imation because it leaves out many-body correlations.^^-) 

In summary, we have found the uniform Hartree-Fock 
ground state of a 2DEG in the presence of spin-orbit 
Rashba coupling, characteristic of asymmetric semicon- 
ductor quantum wells. We present a phase diagram 
where two competing quantum phases are identified, one 
of which shares the Rashba single-particle orbitals, with 
their IP spin quantization axis, and another one which 
has a finite component of the spin in the perpendicular 
direction. This phase possesses a partial spin polariza- 
tion and exists within a window of densities, suggesting 
that a combination of exchange and Rashba spin-orbit in- 
teraction may qualitatively explain experimental results 
obtained in GaAs-AlGaAs heterostructures at low elec- 
tron density. 
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